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Recent measurements of the ^(lS), the ground state of the bb system, show the splitting between 
it and the T(1S) to be 69.5±3.2 MeV, larger than lattice QCD and potential model predictions, 
including recent calculations published by us. The models seem unable to incorporate such a large 
hyperfine splitting within the context of a consistent description of the energy spectrum and decays. 
We investigate whether a perturbative treatment of our potential model (described below) can lead 
, to such a consistent description, including the measured hyperfine splitting, by not softening the 

04 ■ delta function terms in the hyperfine potential. With this modification, we calculate the IS hyperfine 

splitting to be 67.5 MeV, with little effect on the overall fit of our model results. We also present 
, predictions for the 2S and 3S hyperfine splittings. 
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INTRODUCTION 



Recent measurements by the BABAR Collaboration [T|, Q have located the ground state of the bb system, the 
r]b(lS), at a mass of 9390. 8±3. 2 MeV. This value has recently been confirmed by the CLEO Collaboration^. 
Thus, the hyperfine splitting between the IS states is 69.5±3.2 MeV. This splitting is surprisingly large when 
compared with the corresponding splitting in charmonium, given the mass dependence of the leading order 
contribution. In addition, this hyperfine splitting is larger than predictions in recent lattice QCD calculations, 
Q\ ' which range from 52.5±1.5 MeV to 65.8±4.6 MeV [J-tZl' as wen as th° se predicted by recent potential models 
lO ■ (42±13 MeV) A comparison of the data with various model results is presented in the concluding section of 
■ this paper. 

This large hyperfine splitting has recently been investigated as a possible indicator of new physics as a 
means of adjusting the value of the strong coupling parameter and to fix the value of as in an investigation 
other quarkonium states [Toj ] . While further data should clarify these points, in this paper we investigate whether 
£^ | the model we employed in our 2007 paper [ill ], is robust enough to accommodate this new data. 

The results in [Til ] , while providing a good overall fit to the bb spectrum, do not yield such a large hyperfine 
splitting. In that paper, we showed that the perturbative treatment of a model consisting of a relativistic kinetic 
energy term, a linear confining term including its scalar relativistic corrections and the complete perturbative 



one-loop quantum chromodynamic short distance potential was able to reproduce the overall spectrum of the 
bb system as well as its radiative decays, with good accuracy. However, in our model results, the IS hyperfine 
splitting, at 47 MeV, was considerably smaller than the BABAR measurements indicate. 

In our previous calculation, we followed the standard practice of softening the delta function terms in the 
potential [l2|]. This is essential when the entire interaction is treated non-perturbatively in order to avoid 
instability in the numerical calculations. However, in a perturbative treatment, we can perform the calculations 
with the delta function terms unsoftcned and still retain the overall goodness of our fit for the spectrum and 
leptonic decays, while reproducing the correct IS hyperfine splitting, with only minor changes in the potential 
parameters. 
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2. HEAVY QUARKONIUM HYPERFINE POTENTIAL 

The one-loop hyperfine terms arising in the QCD potential are [IH [3] 
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where as is the strong coupling constant in the GR scheme |15l |. 

A variety of approaches to softening the delta function terms have been used (l6j . In most instances, this 
softening is done because the approach to integrating the Schrodinger equation precludes the inclusion of delta 
function terms. In our previous calculation we chose to adopt the quasistatic approximation and softened 
the delta function as 



(2) 



In the interest of clarity, the reevaluation of the hyperfine splittings presented here begins with a variational 
calculation of the wave functions and energies of an unperturbed Hamiltonian. The Hamiltonian consists of the 
relativistic kinetic energy, a linear confining potential and a short-ranged Coulomb-like potential that includes 
the one-loop correction to the strong coupling parameter as ■ This procedure results in an orthogonal set of radial 
wave functions for every orbital angular momentum. These radial functions and the associated energies describe 
the spin-averaged upsilon spectrum. We include the contributions of all v 2 /c 2 and one-loop QCD corrections 
to the the bb potential energy using perturbation theory, retaining any delta function terms that arise in the 
derivation of these corrections. Further details can be found in Rcf. 111. 



The retention of the delta function terms in our calculation is motivated by their use in familiar QED ap- 
plications. The inclusion of delta function contributions in perturbative calculations is necessary to reconcile 
the hydrogen fine-structure results with the Dirac equation, to derive the hydrogen hyperfine splitting, and to 
understand the positronium and muonuim spectra [l8j . Since the wave functions in Ref. (Tl| are derived vari- 
ational^ using a non-singular Hamiltonian, there are no technical difficulties in including delta function terms 
perturbatively. Given the numerous approaches to the calculation of the hyperfine splittings in quarkonia, in- 
cluding the delta function approach of Ref. [To[ , our perturbative evaluation of the hyperfine intervals is relevant 
to the overall discussion. 



RESULTS 



As discussed above, we have recalculated the bb spectrum by retaining the delta function terms in our_ pertur- 
bative calculation. The resulting fitted parameters are shown in Table HI along with those from Rcf. 
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Softened [111 


Unsoftened 


A (GeV 2 ) 


n 1 77+0.006 
U - L ' 1 -0.002 


0.175 


as 


296 +a004 
u.z»u_ 007 


0.295 


m, (GeV) 




5.33 


/i (GeV) 


4.74 


4.82 


fv 


0.00 


0.00 



TABLE I: Fitted parameters for the softened and unsoftened potentials 

As can be seen from Table U the retention of the delta function terms leads to a very minor adjustment of 
the fitted parameters. Our results for the bb spectrum are shown in Table HU It can been seen that the only 
significant changes are in the s-state hyperfine splittings. We calculate these splittings to be: 67.5 MeV for IS, 
35.9 MeV for 2S, 30.3 MeV for 3S. The latter two values are, of course, predictions. 
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m bb (MeVJ 


Softened 


Unsoftened 




Vb(lS) 


9413.70 


9392.91 


9390.8 ± 3.2 


T(15)* 


9460.69 


9460.38 


9460.30 ± 0.26 


xbo(ipy 


9861.12 


9861.39 


9859.44 ± 0.52 


xbi(ipy 


9891.33 


9891.33 


9892.78 ± 0.40 


Xb2{lP)* 


9911.79 


9910.63 


9912.21 ± 0.40 


h b (lP) 


9899.99 


9899.93 




Vb(2S) 


9998.69 


9987.42 




T(2S)* 


10022.5 


10023.3 


10023.26 ± 0.31 


T(1D) 


10149.5 


10149.8 




1 S D 2 


10157.1 


10157.3 


10161.1 ± 1.7 


I s D a 


10162.9 


10163.1 




t 1 T-! 

rD 2 


10158.4 


10158.6 




xbo(2py 


10230.5 


10230.5 


10232.5 ± 0.6 


Xti(2P)* 


10255.0 


10254.8 


10255.46 ± 0.55 


X6 2(2P)* 


10271.5 


10271.2 


10268.65 ± 0.55 


h b (2P) 


10262.0 


10261.8 






10353.0 


10353.1 




l 3 fs 


10355.8 


10355.8 




1 3 F 4 


10357.5 


10357.5 




1 F 3 


10355.9 


10356.0 




Vb (35) 


10344.8 


10333.9 




T(3S') 


10363.6 


10364.2 


10355.2 ± 0.5 


T(2D) 


10443.1 


10443.0 




2 3 D 2 


10450.3 


10450.1 




2 3 Z> 3 


10455.9 


10455.7 




2 1 D 2 


10451.6 


10451.4 




2 3 F 2 


10610.0 


10609.6 




2 3 F 3 


10613.0 


10612.5 




2 3 f4 




1UU14. 




2 : F 3 


10613.2 


10612.7 




^(45) 


10622.8 


10609.4 




T(45) 


10643.0 


10636.4 


10579.4 ± 1.2 



TABLE II: Results for the bb spectrum using softened and unsoftened potentials are shown. Our perturbative fits use the 
indicated states. The value of the 7]b{lS) mass is taken from || and all other data is taken from 

We have also examined the leptonic widths as shown in in Table ITnl Wc find that there is a noticeable increase 
only for the 3S and 4S states. However, the modified results are still compatible with the experiments. 



4. CONCLUSION 

We conclude with a comparison of various modelling approaches with the experimental data. In Table IIV| wc 
show a comparison of the data with the results of representative model calculations: our new (unsoftened) and 
old (softened) models, several recent lattice QCD results, next-to-leading logarithmic perturbative QCD, and a 
QCD-inspired phenomenological treatment. This comparison shows the range of predictions of various modelling 



4 



1 ee (ke V j 


Softened 


Unsoftened 


rjXpt 


TCI "I s ! 


l.OO 


l.OO 


1.04:U ZU U.UIO 


T(25) 


0.61 


0.62 


0.612 ±0.011 


T(35) 


0.46 


0.48 


0.443 ± 0.008 


T(45) 


0.35 


0.40 


0.272 ± 0.029 



TABLE III: The leptonic widths of the T(nS) states using softened and unsoftened potentials, 
approaches. Additional data should provide clarification. 



AHF (MeV) 


Our 


Ref. [11] 


Ref. [4] 


Ref. [5J 


Ref. [6] 


Ref. [7] 


Ref. [20] 


Ref. [10] 


Expt [2] 


AHF(1S) 


67.5±0.7 


47.0 


52.5±1.5 


54.0±12.4 


61±14 


65.8±4.6 


39.5±8 


63.4 


69.5± 3.2 


AHF(2S) 


35.9±0.3 


23.8 






30±19 






35.0 




AHF(3S) 


30.3±0.2 


18.8 












27.6 





TABLE IV: Comparison of s-state hyperfine splittings in various models with experiment. The errors on the hyperfine 
splittings in column 'Our' were obtained by using the Gaussian errors on the parameters as, and /j, associated with 
our fit to the upsilon spectrum. 

We have shown that by not softening the delta function terms which arise in the one-loop quark-antiquark 
potential, we are able to reproduce the surprisingly large hyperfine splitting in the IS level of the bb system. This 
is in contrast to other modeling results, including ours, which soften the delta function terms. In addition, the 
improvement in the hyperfine splittings does not affect either the overall fit of our perturbative calculation to 
the data, nor to our ability to reproduce the leptonic decay widths. We also predict somewhat larger hyperfine 
splittings for the other s-states. 
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